General formulation of transverse hydrodynamics * 
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General formulation of hydrodynamics describing transversally thermalized matter created at the 
early stages of ultra-relativistic heavy-ion collisions is presented. Similarities and differences with 
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the conservation laws as well as the thermodynamic consistency of two-dimensional thermodynamic 
variables characterizing transversally thermalized matter is emphasized. 
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I. INTRODUCTION 

We present a general formulation of the hydrodynam- 
ical model describing transversally thermalized matter 
possibly created at the early stages of ultra-relativistic 
heavy- ion collisions. The concept of purely transversally 
thermalized matter obeying hydrodynamic evolution was 
initially formulated by Heinz and Wong [J 0]. A new 
and different implementation of this idea has been re- 
cently given in Refs. [3, 01 1 where, in contrast to the 
original formulation, it has been shown that this con- 
cept can be consistent with the RHIC data describing 
transverse-momentum spectra of pions and their elliptic 
flow coefflcient V2- 

In this paper we extend the formal results of Refs. 
0, 0] and introduce a general formalism for treating the 
transversally thermalized matter. Similarly to the stan- 
dard relativistic hydrodynamics ^ , our formalism is based 
on a specific form of the conserved energy-momentum 
tensor. We do not restrict our consideration to the case 
of massless particles and the classical (Boltzmann) statis- 
tics, as it has been assumed for simplicity in [1, Q- With 
the given form of the energy-momentum tensor we derive 
the hydrodynamic equations using the thermodynamic 
identities only. We also show how the discussed struc- 
ture of the energy-momentum tensor and other currents 
follow from the explicitly covariant form of the underly- 
ing phase-space distribution function. 

Similarly to Refs. [1, 0] we view our system as a su- 
perposition of non interacting (transverse) clusters. The 
clusters are formed by particles moving with the same 
values of the rapidity. The correlation between the ra- 
pidity y and the spacetime rapidity rj of the particles is 
built in the model by the condition 



y ^ V- 



(1) 
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^ By standard hydrodynamics we always mean the three- 
dimensionally thermalized hydrodynamics of perfect fluid. 



The clusters form two-dimensional (2D) objects, whose 
thermodynamic properties are described by appropriate 
2D thermodynamic variables. One of the advantages of 
our approach, in comparison with the original formula- 
tion by Heinz and Wong, is that our approach is ther- 
modynamically consistent. The 2D thermodynamic vari- 
ables are derived from the formula for the potential 
(valid for massive bosons or fermions) and satisfy stan- 
dard thermodynamic identities. 

The use of the thermodynamic identities together with 
the conservation laws for the energy-momentum tensor 
leads to the entropy conservation. This behaviour is ex- 
pected since the longitudinal free-streaming does not pro- 
duce the entropy and the transverse expansion is entropy 
conserving. As a consequence, our system may be viewed 
as a superposition of 2D perfect fluids. 

Clearly, the microscopic collisions between the parti- 
cles lead to the three-dimensional (3D) isotropization of 
the system, hence the proposed picture of purely trans- 
verse hydrodynamics may be valid only at the very early 
stages of the evolution of matter. Recently, the possible 
2D 3D transition has been addressed in the frame- 
work of the dissipative hydrodynamics in Ref. [Sj] where 
the initial asymmetric pressure relaxes towards the equi- 
librium isotropic form, see also Ref. 

Below we discuss how our formalism may be supple- 
mented with the 2D — > 3D transition, e.g., by apply- 
ing the Landau matching conditions between the purely 
transverse hydrodynamic expansion and standard 3D hy- 
drodynamics. At the transition point the Landau con- 
ditions require local energy-momentum conservation and 
entropy production, which resembles the situation known 
from the dissipative hydrodynamics 

We expect that the early dynamics in relativistic 
heavy-ion collisions is dominated by gluons, whose num- 
ber is not necessarily conserved, hence we restrict the for- 
mulation of the final equations to the case of zero chem- 
ical potential. 

Our main results are presented in Sect. HT] The discus- 
sion of the two-dimensional thermodynamics is presented 
in Sect. IIIII The structure of the phase-space distribu- 
tion function and its moments are presented in Sect. IIVI 
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Below we use the following definitions for rapidity and 
spacetime rapidity, 



1, E+pn 1 t + z 

y = — m — , 77 = — m ■ 



2 E-p\\- 



2 t- z' 



(2) 



which come from the standard parameterization of the 
four-momentum and spacetime coordinate of a particle, 



= (i?,p^,p||) = (mj^ coshy,pj^, mj^ sinhy) . 
— (i, a;^, z) = (r cosh?], a;^, T sinh?]) . 

In Eq. ([3|) the quantity is the transverse mass 



= i/m^ +p2 +p2 



and r is the proper time 



(3) 



(4) 



(5) 



Throughout the paper we use the natural units where 
c—\ and h = 1. 



II. HYDRODYNAMICAL EQUATIONS FOR 
TRANSVERSALLY THERMALIZED MATTER 

A. Energy-momentum conservation laws 

Our main resvilt is that the hydrodynamic equations for 
transversally thermalized matter follow from the energy- 
momentum conservation laws 

d^.T^'' = 0, (6) 
with the energy-momentum tensor of the form 
TM"' ^ [(£2 + uf^U" - P2 {g^"" + VV) ] . (7) 

T 

Here ng is the density of transverse clusters in rapidity 
and r is the proper time defined by Eq. ([5]). 

The quantities £2 and P2 are the 2D energy density 
and 2D pressure. They are scalar quantities determined 
in the local rest-frame of the fluid. The subscript 2 de- 
notes that those quantities do not refer to the 3D volume 
but to 2D transverse area of the cluster. Below, we also 
introduce the 2D entropy density S2, the 2D particle den- 
sity 712, as well as the temperature T and the chemical 
potential /i. In Sect. Illll we show in more detail how 
these quantities are introduced. The important point is 
that all 2D thermodynamic variables are defined in the 
consistent way and satisfy well-known thermodynamic 
identities. For example, for vanishing chemical potential 
we have 



and 



£2 + P2 = TS2 



de2 = Tds2, dP2 — S2dT. 



(8) 



(9) 



In Eq. ([7]) the four vectors and are defined by the 
equations: 



= (mq cosh ry, ita;, ttj,, 77,0 sinh?7). 



(10) 



(11) 



(sinhry, 0, 0, cosh 77), 
and 77^,^2;, Uy are the components of the four vector 

= {u",u^,0) ^ {u",u,,uy,0) , (12) 



which describes the four-velocity of the fluid element in 
the rest frame of the cluster to which this element be- 
longs. The four-velocity u'^ is normalized to unity 



7^77^ ' 



77^ = 1. 



(13) 



The four- vector describes the four- velocity of the fluid 
element and plays a standard role of the hydrodynamic 
flow. It may be obtained from 77^ by the Lorentz boost 
along the z axis with rapidity rj (in other words, com- 
bines the motion of the fluid element in a cluster with 
the motion of the cluster). The appearance of the four- 
vector is a new characteristic feature of transverse 
hydrodynamics related with the special role of the longi- 
tudinal direction - in the rest frame of the cluster we have 
y = (0, 0, 0, 1). The four vectors W and V satisfy the 
following normalization conditions: 



W'U^ = 1, 



y^K = -1, c/^K = 0. 



(14) 
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The presence of the term V^V^ in ([7]) is responsible for 
vanishing of the longitudinal pressure. In the local rest- 
frame of the fluid element, where we have = (1, 0, 0, 0) 
and Vf" = (0, 0, 0, 1) one flnds 



(15) 



The fact that T^^ = indicates that there is no interac- 
tion between the clusters. Certainly, such situation can- 
not last long and we expect that a transition from trans- 
verse to standard 3D hydrodynamics must take place 
within at most few fermis after the collision. 



B. Entropy conservation 

By straightforward calculation, we can show that the 
energy-momentum conservation laws ^ lead to the en- 
tropy conservation, namely, we flnd 

U.d^T^" = Td^S^ = 0, (16) 

where the entropy current is dcflned by the expression 

(17) 



5^ = ^s2C/^ 
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In Sect. |TV]we show that both (O and ([TTI) follow from 
the same phase-space distribution function, which real- 
izes our main assumption of transverse thermalization 
and longitudinal free-streaming. Thus, the definition of 
the entropy current is consistent with the definition of 
the energy- momentum tensor 

In the derivation of Eq. (fTB|) we used the relations 



"0, 

■"0, 



(18) 
(19) 



which follow directly from the definitions ([5]), (fTU]) and 
PT|) . It is worth emphasizing that Eq. has the 

structure known already from the standard relativistic 
hydrodynamics. In both cases the entropy conservation 
(the adiabaticity of the flow) is a direct consequence of 
the energy-momentum conservation. The situation is dif- 
ferent in the non-relativistic hydrodynamics where the 
entropy conservation typically appears as an additional 
condition supplementing the standard Euler equation. 



C. Euler equations 

The use of the entropy conservation ([TB|) in ^ leads 
to the equation which may be treated as the analog of 
the Euler equation in our case. 



(20) 



By performing projections of Eq. (|20p on the four- vectors 
C/y and V^, which yield identically zero, we verify that 
only two out of four equations in pop are independent. 
In this calculation we use the relations 



0, 
0, 



(21) 
(22) 
(23) 



which, similarly as Eqs. ^TE\i and (|19p . follow from the 
definitions ([5)1. (fTU)) and (fTT|) . 

From the physical point of view the fact that pp]) con- 
tains only two independent equations is expected. Since 
the longitudinal motion is fixed, the Euler equation may 
specify only the dynamics of the two transverse compo- 
nents of the velocity, i.e., it determines only the trans- 
verse expansion of matter in a cluster. 

Altogether we have at our disposal three equations, 
one in (|16p and two in (|20p . for four unknown functions: 
S2, T, Ux and Uy. To close this system we have to add 
the equation of state, which specifies the temperature 
dependence of the entropy. This situation reminds the 
case of the standard relativistic hydrodynamics where we 
have four equations for five unknown functions and also 
the equation of state must be added to close the system 
of equations ^ . 



D. Cylindrical coordinates 

The three independent equations of transverse hydro- 
dynamics may be conveniently written in the cylindrical 
coordinates. In this case we introduce the distance r and 
the azimuthal angle <p defined by the equations, 



= tan ^{ry/r^), 



(24) 



where x±_ = [rx.ry) 
the fluid velocity. 



In the similar way we parameterize 



V cos(a -I- 1 



V sin(Q! -|- I 



UqVx 

uoVy 



u± cos(a -|- cj)), 
u± sin(a -I- (j)). 



(25) 
(26) 
(27) 



Here v is the transverse flow, uq — (l ~ v"^) ^ ,u± = uqv, 
and the dynamical angle a describes deviations of the 
flow direction from the radial direction. 

With the help of such variables the hydrodynamic 
equations may be written explicitly in the following form 

d d d 

— {rs2Uo) + — {rs2UQVC0sa) + — (s2Uo«sina) = 0, 
OT or a<p 

d d d 

— (rTuov) + r cos a— (Tuq) + sm a— (Tuq) = 0, 



1 UnV I — 



V sma 



sma 



dT cos a dT 



dr 



d(j) 



= 0. (28) 



In the derivative d/dr denotes the total derivative 
with respect to time 



d 

d^ 



d_ 

dT 



V cos a 



d V sin a d 



dr 



(29) 



We note that Eqs. pgj) have exactly the same struc- 
ture as Eqs. (20) in Ref. The important point is, 
however, that they were derived here directly from the 
energy-momentum tensor ([SJ with the help of the ther- 
modynamic identities only. No specific dependence of the 
thermodynamic variables on the temperature was used. 



E. Breaking of boost-invariance 

So far we treated the cluster density no as a constant 
parameter. One may check, however, that rip may de- 
pend on the spacetime rapidity rj. To see this point, we 
write the the energy-momentum tensor as the product 



(!)■ 



(30) 



^ This is so in the case of vanishing baryon potential. For non- 



zero baryon potential we have one more thermodynamic variable 
(baryon number density) and one more equation (the baryon 
number conservation law). 
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The energy-momentum tensor T^,'^^ agrees with our pre- 
vious form for the case tiq = 1 and is conserved, 



8 T'"' — 



(31) 



It turns out that the conservation of the tensor T^'^'j' im- 
pUes the conservation of the tensor T^". To prove this 
point we consider the expression d^T^''' . Since the tensor 
T^^ is conserved, we find 



Tl'^'^d.noiv). (32) 



With the cxphcit form of Tr^-. we get 



T 



(33) 



Since the operator U^df^ does not contain the derivative 
with respect to ?/, the non-zero contributions on the right- 
hand side of Eq. (|33p may come only from the terms 
in the second hne with v = ot v = 3. The exphcit 
calculation shows, however, that these terms yield zero, 
hence the tensor T^'' is conserved. 

The fact that no may depend on rj indicates that our 
formalism is not necessarily boost-invariant. Actually, 
a more detailed analysis shows that our three hydrody- 
namic equations describe the expansion of matter in a 
single cluster and different initial conditions may be ap- 
plied for different values of rj. Such situation is of course a 
direct consequence of our main assumption that clusters 
are independent. 



Landau matching conditions at 2D 
transition 



3D 



The 2D — !■ 3D transition may be described by assuming 
the Landau matching conditions at the transition point 



(34) 



where is the standard energy-momentum tensor of 
relativistic hydrodynamics of perfect fluid 



T^^ = {e + P)U^'U'' - Pg^'r (35) 
Substitution of Eqs. ^ and ([23) into Eq. ([M]) gives 



T 



(36) 



This condition specifies nothing else but the local conser- 
vation of energy and momentum at the transition point. 
In addition, Eq. p6p must be supplemented by the con- 
dition that the entropy increases 



no 



S2 < s. 



(37) 



Here s is the 3D entropy density corresponding to the 3D 
energy density e. 

We note that a similar procedure of Landau matching 
has been recently used in Ref. [3], where the 3D free- 
streaming in the early stage is followed by the 3D boost- 
invariant hydrodynamics. 



III. THERMODYNAMICS OF 
TWO-DIMENSIONAL SYSTEMS 

The starting point of our thermodynamic considera- 
tions is the formula for the potential Q of non-interacting 
bosons (upper signs) or fermions (lower signs). In the 
case of two-dimensional systems it has the following form 



n{T,V2, ti) ^ ±iygTV2 



(27r)2 



In 1 T e 



(M-m_L)/T 



(38) 

The number of particles is not conserved in our approach 
and we assume /i = (however, this can be done only in 
the end of the thermodynamic calculations). Since the 
motion of particles is confined to a plane, their energy is 
equal to their transverse mass defined by Eq. (|4]). 

The quantity V2 is the area of the plane where the par- 
ticles move and Vg denotes the internal degrees of free- 



dom. For gluon dominated system i>g — 16. Eq. 
allows us to define the particle density n2 = N2/V2, pres- 
sure P2, entropy density §2 — S2/V2, and energy density 
£2 = i?2/V2. Those quantities follow from the well-known 
thermodynamic relations: 



A^2 - 
P2 = 
S2 = 



dV2/Ta 

on 



V2' 



and 



E2+P2V2^TS2+tiN2. 



(39) 
(40) 
(41) 

(42) 



In our calculations the last equation is most often used 
in the form valid for the densities, see Eq. ([8]). In the 
case /i = Eqs. ^ - ^ yield: 



iV2 
P2 
S2 

E2 



(27r)2 ^' 
^ I (2^2^ 



g, 



(43) 
(44) 



= -l^qV2 



^9V2 



(2^)2 



[<?ln5T(l±5)ln(l±5)], (45) 



fpi 
i2ny 



■m±g, 



(46) 
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where we have introduced the equihbrium distribution 
function 

Here again the upper (lower) signs refer to bosons 
(fermions). 



IV. LORENTZ STRUCTURE OF THE 
PHASE-SPACE DISTRIBUTION FUNCTION 
AND ITS MOMENTS 

In this Section we analyze in more detail the Lorentz 
structure of the phase-space distribution function which 
leads to the form of the energy-momentum tensor dis- 
cussed in Sect. II. 



A. Massless fermions and bosons 

In the case of massles fermions or bosons, the integrals 
and are analytic and yield: 

712 — — , £2 — — ^_ — (fermions) (48) 



24 



47r 



and 



n2 = ^— , £2 ^ ' (bosons). (49) 

In both cases we also find 



1 2 1 



(50) 



The entropy density may be obtained from Eq. ([5]) with 
the appropriate substitutions. 



Distribution function 



According to the original formulation of the model in 
Ref. 0], we assume that the phase-space distribution 
function F may be factorized into the longitudinal and 
transverse part 



(57) 



In the global equilibrium, the function g agrees with the 
equilibrium distribution (|47p . In the local equilibrium, 
assumed in the case of the transverse hydrodynamic ex- 
pansion, it has the form 



(58) 



where 



B. Classical limit 

In the case of the classical (Boltzmann) statistics the 
analytic formulas for the densities of thermodynamic 
quantities may be easily obtained also for finite masses, 



n-2. 

£2 



^(m + r)e-™/^, 
27r 



3mT + 3T2]e-™/'^, 



2tt' 

^[r2 + (m + T)2]e-WT^ 
27r 



(51) 

(52) 
(53) 
(54) 



In the limit to ^ these results agree with the expres- 
sion first given in Ref. 0]. we note that in the case of 
the Boltzmann statistics, the equation of state is of the 
well-known form, P2V'2 = N2T. With the help of the 
thermodynamic identities de = Tds2 and dP2 = S2dT we 
find the formula for the sound velocity, 



dP2 
de2 



S2dT dlnT 



Tds2 dlns2 
T{w? + 3mT + 3T^) 
+ 3m^T + 6toT2 + 6T3 ' 



(55) 



In the limiting cases one obtains the expected results. 



lim = -, 

m^o " 2 



lim ci = —. 

T^oo 2 



(56) 



p'^Uf^ = m_LMo cosh(7; - rj) - ■ Uj_. 



(59) 



The longitudinal part of the distribution function is given 
by the expression which implements the condition ([T|). 



S{y - rn) no 

f\\ = '^o = —S (p^V^) 

TO_Lr T 



(60) 



Combining Eqs. ([57|) . (l58|) and ([60|) we obtain the form 
of the phase-space distribution function which explicitly 
emphasizes its scalar character. 



F=!^S{p-V)9{p-U). 

T 



Particle current 



(61) 



The particle current iV^ is defined as the first moment 
of the distribution function, 



^J'^P'SiP-V)9iP-U). (62) 



The Lorentz structure of the distribution function implies 
that iV'' may be written as the combination of the four- 
vectors and J7^, namely 

m ^ + bU''. (63) 

The coefRcients a and b are obtained by the projections. 



-a, U^Nf' = h, 



(64) 
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and, as scalar quantities, they may be calculated 
in the local rest frame of the fluid element where 
L/^ = (1,0,0,0) and V = (0,0,0,1). The explicit cal- 
culations give 



^^V(2nf J -^^'5(p||).g(i?) - — n^. (66) 

Here we used the definition of the two-dimensional par- 
ticle density following from Eq. (|43|) . Combining all the 
results of this Section we find the general structure of the 
particle current, 

iV^ = ^naC/'^ (67) 



C. Energy-momentum tensor 



where we used p^p^ = m? ■ In the local rest frame we 
find 

n^Vg f (Pp^ 
T J (2tt)^ m± 

= !^(£2-2P2). (72) 

T 

Here we used Eqs. (01]) and (US]). We note that for 
massless particles the trace vanishes, T^^ — 0, and we 
find £2 = 2P2- This in turn gives — 1/2. In the next 
step we calculate the projection of the energy-momentum 
tensor on the four- velocity If^, 

= I^2t^Yt J P^'^^9('m±) = —£2- (73) 
Eas. (|70| - ((73|) lead us directly to our basic definition ((T]). 



The energy-momentum tensor is the second moment 
of the distribution function. It should have the following 
Lorentz structure, 

T^. ^ a'U^U,+b'g^,+c'V^V, + ^iU^V, + U,V^). (68) 

Here, besides the symmetric combinations of the four- 
vectors and y , the term proportional to the metric 
tensor is also required. Simple algebra gives 

T^^ = a' + 4b'-c', 



T^'^Uf^U^ = a' + b', 
T'^'Vf.V^ = -b' + c 



d' 



(69) 



The coefficients in the decomposition (|68p may be again 
calculated in the local rest frame of the fluid element. 
This procedure yields immediately the two constraints: 
d' ~ and c' = b'. Thus, Eqs. (p^ are reduced to the 
form 



T^^ = a' + 36', 
T^-'U^U, = a' + b', 
c' = b' 
d' = 0. 

The trace of the energy-momentum tensor equals 
noVg f d^p 



(70) 



(27r 



'f-J'^mH{p-V)g{p-Ul (71) 



D. Entropy current 



Finally, we consider the entropy current. 



'(27r)2r 



(74) 



X J S{p ■ V)g{p ■ U)[ Hg{p • f/)] - 1 ]. 

The meaning of this definition is that we use the Boltz- 
mann expression to calculate the entropy of a single clus- 
ter and then the sum over clusters is performed. Similarly 
to the case of the particle current we express the entropy 
current as a linear combination of the four-vectors If^ 
and 

S^" = a"V>' + b"U>'. (75) 
In the analogous way to the cases studied above we find 



a — V, b — — S2 

T 

hence the definition p?]) is recovered. 



(76) 



V. CONCLUSIONS 

In this paper we have introduced the general formalism 
for hydrodynamical description of transversally thermal- 
ized matter possibly created in the early stages of ultra- 
relativistic heavy-ion collisions. We have showed how the 
hydrodynamic equations follow directly from the general 
form of the energy-momentum tensor. We discussed the 
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entropy conservation law and the structure of the un- 
derlying phase-space distribution. The thermodynam- 
ics of two-dimensional systems has been analyzed and 
the generally covariant structure of the model has been 
emphasized. After first successful attempts to describe 
physical observables studied at RHIC [3| we hope that 



the presented formalism will find further applications to 
describe early stages of relativistic heavy-ion collisions. 
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